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Supersymmetric theories with gravitino dark matter generally do not allow the high reheating 
temperature required by thermal leptogenesis without running afoul of relic abundance or big bang 
nucleosynthesis constraints. We report on a successful search for parameter space that does satisfy 
these requirements. The main implication is the near degeneracy of the gluino with the other 
neutralinos in the spectrum. The leading discovery channel at the LHC for this scenario is through 
mono jet plus missing energy events. 



GRAVITINOS FROM HIGH REHEATING 
TEMPERATURE 



Our goal is to investigate the supersymmetry spectrum 
that allows gravitinos to be the lightest supersymmetric 
partner (LSP) dark matter of the universe and also allows 
thermal leptogenesis to explain the baryon asymmetry, 
all while retaining the successful predictions of big bang 
nucleosynthesis. We explain each of these in turn, high- 
lighting the potential sources of conflict between them, 
and finally settling upon an explanation requiring a neu- 
tralino next-to-lightest supersymmetric partner (NLSP) 
degenerate with the gluino, and then investigating its 
consequences for the Large Hadron Collider (LHC). 

To begin we note that in gauge mediated supersym- 
metry [T] the lightest supersymmetric particle (LSP) is 
generically the gravitino, with mass ranging from the eV 
scale to tens of GeV. There are many nice features of 
this model that we do not detail here, but suffice it to 
say that it is a powerful and viable organizing principle 
for the superpartner spectrum, and motivates our inter- 
est in the gravitino as the LSP. The identity of the next 
to lightest supersymmetric particle and the details of the 
spectrum are model dependent (see e.g., [5]). Being the 
LSP and stable, the gravitino is the leading candidate for 
dark matter in these models, apart from the possibility 
of very long-lived particles in the messenger sector [3J. 
However, a light gravitino is at most Warm Dark Matter 
and not favoured by structure formation; indeed, present 
observations already constrain its mass to be above a few 
keV if it decouples as a relativistic thermal relic [3] . 

This is not the only source of gravitinos in the early 
Universe, as they do not have to reach thermal equilib- 
rium densities to be cosmologically important. On one 
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hand, scatterings in the thermal plasma produce grav- 
itinos with the abundance proportional to the reheating 
temperature after inflation [SHI]: 
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where M r denote physical gaugino masses and the coef- 
ficients 7 r depend on the ratios of the gauge couplings at 
the reheating scale and the scale of the physical gaugino 
masses: with the 1-loop RGE for the gaugino masses, 
the values of 7 r can be evaluated for Tr = 10 9 (10 7 ) GeV 
as 73 = 0.48 - 0.56(0.62 - 0.74), 72 = 0.57(0.54), 
7x = 0.22(0.17), where the range for 73 corresponds 
to the gluino masses ranging from 200 to 900 GeV 
[8J. We have only included the production of the gold- 
stino component of the gravitino, which dominates for 
«tNLSp/™3/2 > O(10). On the other hand, gravitinos are 
also produced in the gravitational decays of the NLSP, 
but for ^nlsp^ 2 « 1 or TI3 12 /^nlsp "C 1 these de- 
cays are a negligible source of gravitino dark matter; 
moreover, too high a fraction of such a nonthermal and 
warmish dark matter component can cause too much era- 
sure of the cosmic structures at small scales [5]. Other 
contributions to the gravitino abundance can arise from 
inflaton decay [TO] or from the reheating process [TJJ , but 
they are more model dependent and we will not discuss 
them further. 

Thus, the gravitino abundance is largely determined 
directly by the reheating temperature ~ Tr, as sug- 
gested by eq. . If this were the only way the reheating 
temperature affected the scenario, one could contemplate 
a simple explanation for the cold dark matter by tuning 
Tr to achieve the required £Iq. 

However, there are other implications to the choice of 
Tr that must be considered. Thermal leptogenesis re- 
quires a rather high Tr to be successful, which may yield 
too much dark matter unless the gravitino mass is lifted 
to higher values (Qg ~ Tr/to 3 / 2 , assuming tonlsp is 
fixed). However, higher gravitino mass means a slower 
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NLSP decay, which follows a normal thermal relic history 
and then dumps its decay energy into NLSP — > G + X af- 
ter big bang nucleosynthesis (BBN). BBN compatibility 
depends on both the number density of the NLSP tinlsp 
and its decay lifetime tnlsp an d other quantities (such as 
the decay branching fractions, the gravitino mass, etc.). 
Thus, there are strains and potential incompatibilities 
when requiring compatibility among gravitino dark mat- 
ter, thermal leptogenesis and big bang nucleosynthesis. 

In the next two subsections we shall describe in more 
detail the constraints that arise from requiring both suc- 
cessful thermal leptogenesis and compatibility with BBN. 
We shall then explain in sec. |I C| our emphasis on sce- 
narios with neutralino NLSP. In sec. [TT] we determine 
the maximum allowed reheating temperature, which is 
wanted by thermal leptogenesis, consistent with all the 
constraints. We show how this value depends on other 
parameters of the theory. In sec. |III| we discuss the impli- 
cations of the resulting parameter space for the Tevatron 
and Large Hadron Collider. We summarize our conclu- 
sions in sec. II VI and make some additional final remarks. 



A. Thermal Leptogenesis Requirements 



on a light gravitino LSP being the dark matter while 
being consistent with BBN, even if one tries to make it 
as large as possible by requiring that the masses of the 
gauginos and the NLSP are not too far apart [5]. That 
consistency is the subject of the next subsection. 

Before closing this discussion, we remark on a few 
caveats to what was said above. It has been argued that 
neglecting quantum effects in the dynamics of leptoge- 
nesis introduces theoretical uncertainties as large as an 
order of magnitude [16] . Furthermore, the bounds dis- 
cussed above do not apply in models with degenerate 
masses of the right handed neutrinos, see e.g., [TT 1 fTH ] . 
or in models with large cancellations in the seesaw mass 
formula |2"U] . However, these solutions are of some- 
what less interest to us here since they involve a degree 
and type of fine tuning that we wish to try to avoid. We 
also note that the reheating temperature may be lowered 
in models of nonthermal leptogenesis, see e.g., [3TJ[22], or 
in soft leptogenesis, see [23H25] , Such models, however, 
require arranging for additional interactions, e.g. a cou- 
pling between the inflaton and a right-handed neutrino 
or a coupling of the Higgs doublet to the leptonic com- 
ponent of the messenger field. But in this approach the 
attractive feature of independence of initial conditions is 
lost - another consequence that we wish to avoid here. 



Generation of the baryon asymmetry through ther- 
mal leptogenesis remains a theoretically attractive and 
experimentally viable possibility, as it only uses parti- 
cles (right-handed neutrinos) and interactions (neutrino 
Yukawa couplings) already present in the seesaw mod- 
els explaining the smallness of the neutrino masses (for 
a review, see |12|). Putting it simply, the lepton asym- 
metries, subsequently converted into baryon asymmetry, 
are produced in the CP violating decays of the lightest 
right-handed neutrinos; it is usually assumed that these 
neutrinos had been previously produced in scatterings in 
the thermal plasma. So successful thermal leptogenesis 
gives a lower bound on the mass of the lightest neutrino, 
M^ m ' ) , in usual seesaw models (with hierarchical masses 
of the right-handed neutrinos), which can be translated 
into a lower bound on the reheating temperature after 
inflation: T R min) « M^f n) /5 (T^ min) « M { ™ n) ) in the 
so-called strong (weak) washout regime. 

Interestingly, the lower bound T R min) =1.9- 10 9 GeV 
[T3] found in the model with the largest possible CP 
asymmetries and optimal (small) washout approximately 
coincides with values of T R obtained in generic models 
predicting the correct baryon asymmetry, found with the 
Monte Carlo Markov chain techniques [2]. However, 
this bound on the reheating temperature is not an abso- 
lute one. If the initial conditions for leptogenesis include 
a thermal distribution of the lightest right-handed neu- 
trino, the minimal reheating temperature is 2.5 • 10 8 GeV 
[15] . This is the uncomfortably high reheating tempera- 
ture with respect to gravitino abundance that we referred 
to above in the introduction. Such a high Tr puts strain 



B. BBN Consistency 

Since the abundance of dark matter (in our case con- 
sisting of gravitinos) is fl^h 2 = 0.110 ± 0.006 [55], sub- 
stituting the minimal reheating temperature J^ min ) dis- 
cussed above into ([lj shows that the gravitino mass 
consistent with the dark matter abundance is at least 
0(1 GeV) in the most optimistic case of the gaugino 
masses degenerate with ttinlsp- With such gravitino 
masses the NLSP lifetime, 



tnlsp = 




for m 3 / 2 /"iNLSP ^ 1) easily exceeds the duration of the 
big bang nucleosynthesis. NLSP decays taking place dur- 
ing or after BBN, can alter its successful predictions if 
the relic abundance of the NLSP and/or the hadronic 
branching fraction in the NLSP decay is large enough 
[271438] . Furthermore, if the NLSP is charged, it can bind 
with nuclei, which facilitates the production of 6 Li [591 - 
41 . This introduces a tension between successful thermal 

leptogenesis, which requires Tr > T R min \ and gravitino 
Dark Matter. For a given MSSM spectrum and known 
f^NLSP^ 2 this tension can be translated with ([I]) and Q 
into a lower bound on tnlsp an d an upper bound on 
m 3 / 2 . 

The BBN bounds are weaker and more easily satis- 
fied for shorter lifetimes, see also [12TI45] . and disappear 
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below 0.1 s. Therefore, there has been an effort to iden- 
tify the NLSP candidates for which the BBN bounds are 
weaker than usual, hence allowing for relatively heavy 
gravitino DM and a high reheating temperature. Several 
solutions have been proposed, involving cither a reduc- 
tion of the NLSP relic abundance compared to the generic 
case, see e.g., [46448] . or the suppressing of the energy 
released in the NLSP decay kinematically, thanks to ex- 
tremely small mass splitting between the NLSP and the 
gravitino LSP [49] . 

The type of NLSP also plays an important role and 
changes the BBN bounds. One of the most studied is 
the stau: it naturally is the NLSP in minimal models 
of gauge mediation with a large messenger number and 
a high scale of supersymmetry breaking (the latter fea- 
ture also predicts a heavy gravitino), and in its decay 
few energetic hadrons dangerous to BBN are produced. 
Nevertheless, the stau is charged and it is constrained 
by bound-state effects, so for a typical stau relic abun- 
dance, a reheating temperature larger than O(10 8 GeV) 
is excluded, even with a compressed spectrum of stau 
and gaugino masses [8]. Stau relic density can also be 
suppressed thanks to a large left-right mixing in the stau 
sector. This effect has been studied in the context of the 
CMSSM and possible reheating temperatures as large as. 
~ 10 9 GeV for fi < gS] and O(10 8 GeV) for /x > [5D] 
were found. There is also a possibility that the stau an- 
nihilation cross section is enhanced by a Higgs pole [5T] . 
All these solutions require a fine tuning among the soft 
supersymmetry breaking mass parameters. 

Sneutrinos as the NLSP easily evade the BBN bounds 
even for a high reheating temperature suitable for ther- 
mal leptogenesis, as long as their masses do not exceed 
200-300 GeV [5H [53] . Hovewer, arranging for a sneu- 
trino NLSP requires a strong degeneracy between the 
soft supersymmetry breaking mass parameters for the 
superpartners of the left- and right-handed leptons (see 
[5] for gauge mediation) or non-universal Higgs masses 
(sec [54 56 for gaugino and gravity mediation) . 

Therefore, our interest turns to a relatively unexplored 
option, that of thermal leptogenesis with a neutralino 
NLSP and a general spectrum of supersymmetric par- 
ticles, as most recent studies on neutralino NLSP were 
done in the context of CMSSM with a gravitino LSP. In 
the next section we shall describe the challenges of this 
option, then do the quantitative work in sec.[ll]to show 
that it can work under some circumstances, and then 
describe the LHC implications in sec. |III[ 



C. Leptogenesis with Neutralino NLSP 



Bh ranges from 3 to 50 percent for neutralino masses be- 
tween 100 and 1000 GeV 07]). Secondly, the bino NLSP 
usually has a large relic abundance and it has been con- 
sidered mostly in the context of models with universal 
gaugino masses, i.e. with M r /g% independent of the gauge 
group index r at the scale of supersymmetry breaking 1 . 
With the gluino mass approximately 5 times larger than 
the bino mass at the electroweak scale, it is clear from 
eq. ([I]) that the resulting reheating temperature is smaller 
by an order of magnitude with respect to the case of lit- 
tle or no hierarchy between the gaugino masses. As will 
become evident in the figures of sec.[TT] the typical maxi- 
mum reheating temperature for the bino NLSP with uni- 
versal gaugino masses reaches only 10 5 — 10 6 GeV, which 
is much too low for successful thermal leptogenesis needs. 

Reconciling models of gauge mediated supersymmetry 
breaking with thermal leptogenesis requires either low- 
ering the leptogenesis temperature or relaxing the BBN 
bounds. As discussed above, this can be achieved sim- 
ply by reducing the relic density of the NLSP e.g. by 
means of coannihilations with more strongly interacting 
particles (which also evades strong cosmological bounds 
on the presence of metastable charged or strongly inter- 
acting particles). Reducing the NLSP relic density by 
coannihilations remains largely unexplored phenomeno- 
logically, expecially in this context, but has been ex- 
plored for the case of bino-stau coannihilation in the 
CMSSM [3T1 133 [Ml 13S1 EE] or in more general mod- 
els g7]. 

In the following, we would like to focus on the most 
promising case of neutralino/gluino mass degeneracy. We 
determine if the reheating temperatures allowed in such 
a scenario are consistent with thermal leptogenesis and 
what fine tunings may be necessary. This part of our dis- 
cussion depends only on the assumption that gravitinos 
are the dark matter particles, and we assume arbitrary 
masses of the supersymmetric particles, without impos- 
ing constraints from specific theoretical models. Note 
that for bino and wino NLSP, the gaugino masses play a 
domininant role in the gravitino and NLSP abundances, 
and a degeneracy helps to reduce both. Therefore the 
degenerate spectrum can be considered the best case sce- 
nario. 



II. MAXIMIZING THE REHEATING 
TEMPERATURE 

Whether a given set of the parameters leading to grav- 
itino LSP and a neutral NLSP is consistent with the 
BBN bounds depends mainly on 4 quantities: the NLSP 



At first sight, neutralinos as candidates for the NLSP 
appear much worse than sleptons or sneutrinos |34j . 
First, hadrons are often found in the decay channels of 
the neutralinos (roughly, the hadronic branching fraction 



At 1 loop, M r /g^. is a renormalization group invariant; for con- 
cretness, we impose the universality condition at the boundary 
scale 10 15 GeV. 
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mass, its relic abundance, its hadronic branching frac- 
tion and its lifetime (which can be traded for 7/13/2 with 
the use of ([2])). Refs. [371 [35] nn d the excluded regions 
on the (tnlsPj ^nlsp^ 2 ) plane for different values of the 
hadronic branching fraction and two values of the NLSP 
masses, 100 and 1000 GeV. Up to a moderate shift in 
the allowed ^nlsp^ 2 > the excluded regions are very sim- 
ilar for both these masses. In order to apply the BBN 
bounds for a general set of parameters, we interpolate 
(linearly on a logarithmic scale) between the results of 
[571 I3"5] , constructing the maximal allowed NLSP abun- 
dance, f^£gpft, 2 (mNLSPj r NLSP: Bh). Using this function, 
for a given MSSM spectrum, we can calculate fiNLSP^ 2 , 
find the maximal 777,3/2 for which Onlsp^ 2 < ^nlsp^ 2 
and, from ([!]), find the maximal allowed reheating tem- 
perature. 

In general, ^nlsp^ 2 decreases with growing Bh, but 
it is a non-monotonic function of tnlsp : it consists of 4 
convex parts, reflecting the bounds for the abundance of 
4 He, D, 6 Li and 3 He. Since for increasing NLSP mass the 
predicted ^nlsp^ 2 an d Bh increase, at certain values of 
777NLSP) we may find the BBN-allowed regions (around 
local maxima of ^nlsp^ 2 as a function of tnlsp) close. 
This makes the bounds for 777,3/2 and Tr discontinuous 
functions of ttinlsp- 

In the case of no coannihilations, we calculate the 
NLSP relic density with the computer code mi cr Omegas 
[571 155) . In the case with neutralino/gluino degener- 
acy we include the coannihilations (taking into account 
nonperturbative contributions) using the prescription de- 
scribed in detail in Appendix A. For hadronic branching 
fractions for bino, wino and higgsino NLSP we use the 
results obtained in |47j , while we take it equal to 1 in the 
gluino NLSP case. 



A. Bino NLSP 

Our results for bino NLSP are shown on Figures [T] 
[2] and |3J The maximal allowed 777,3/2, corresponding in 
practice to the largest reheating temperature, can be di- 
rectly obtained from the BBN bounds. Figure [T] shows 
the predictions for Onlsp^ 2 for three cases: without de- 
generacy with gluino, with 10% and 1% degeneracy, to- 
gether with the BBN bounds for three different values 
of the gravitino mass, 7773/2 = 0.1, 1, 10 GeV and two 
values of Bh — 0.01, 1. In the case of non-degeneracy, 
it is necessary to know the remaining susy spectrum for 
the purposes of computing the bino relic abundance. For 
that, we have chosen the spectrum of minimal gauge me- 
diation with one 5 + 5 messenger pair and the messen- 
ger mass of 10 15 GeV. For such a spectrum the main 
contributions to the bino pair annihilations come from t- 
channel slepton exchanges and the bino-to-slepton mass 
ratio is approximately 0.4. Since the bino annihilation 
cross-section is proportional to l/m~, by increasing this 
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FIG. 1: Gray dashed lines show the BBN bounds on 
{rtig, fiNLSP^ 2 ) pl ane for three different values of the grav- 
itino mass, m 3 / 2 = 0.1, 1, 10 GeV and two values of Bh = 
0.01, 1 Also shown are predictions for SInlsp^ 2 for the non- 
degenerate NLSP case (labeled NO) and with NLSP/gluino 
degeneracy of 10% and 1%; thick black solid (dashed) lines 
correspond to the bino (wino) NLSP. 

ratio (e.g. increasing the numbers of messengers in min- 
imal gauge mediation), one can suppress the resulting 
Qnlsp^ 2 - However, even with 1% bino/slepton mass de- 
generacy, Onlsp^ 2 is smaller only by a factor of ~ 20 
compared to our reference case, so it is still larger than 
what we obtain with 10% bino/gluino mass degeneracy. 

Figure [2] translates these results into the largest al- 
lowed m 3 /2 for each of the three cases and Figure [SJ shows 
the largest allowed reheating temperature for each case. 
Since the reheating temperature depends on the pattern 
of gaugino masses, we assume universal masses for the 
case without degeneracy and a pattern Mi : Mi : M3 w 1 : 
2:1 for the degenerate cases. 

We see that with 1% bino-gluino mass degeneracy, we 
are able to reach a reheating temperature as high as a few 
10 9 GeV with < 300 GeV. Even with only 10% bino- 
gluino mass degeneracy, we are able to reach Tr > 0.7 • 
10 s GeV for bino masses all the way up to 1 TeV. 



B. Wino NLSP 

Our results for wino NLSP are also shown on Fig- 
ures [T] [2] and [3] in the same way as for the bino NLSP 
case. The only difference is instead of the universal 
gaugino mass pattern, which always gives the bino as 
the lightest gaugino, we use the spectrum arising in 
models with anomaly mediated supersymmetry breaking 
(AMSB) [59, 60J. This serves as a good illustration of a 
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FIG. 2: BBN bounds on (mg, 1713/2) plane for the non- 
degenerate NLSP case (labeled no B — g (w — g) degeneracy) 
and with NLSP/gluino degeneracy of 10% and 1%; thick black 
solid (dashed) lines correspond to the bino (wino) NLSP 
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FIG. 3: BBN bounds on (m B ,Tn) plane. Solid lines corre- 
spond to the bino NLSP case with universal gaugino masses 
and to bino/gluino degeneracy of 10% and 1% with Mi : 
Mi : M 3 w 1:2:1 (labeled 1-2-1). Dashed lines corre- 
spond to the wino NLSP case with AMSB gaugino mass spec- 
trum and to wino/gluino degeneracy of 10% and 1% with 
Mi : M 2 : M 3 « 2 : 1 : 1 (labeled 2-1-1). 



theory approach giving the wino as the lightest gaugino. 

Models with nondegenerate wino NLSP allow for a re- 
heating temperature from a few 10 6 GeV for the AMSB 
gaugino spectrum, with m g ~ 10 m^, to a few 10 r GeV 
for milder mass hierarchies (given by the 10% w — g 
degeneracy line with almost the same ^nlsp^ 2 as the 



nondegenerate case). Although such Tr is larger than 
in the models with the generic bino NLSP, a degener- 
acy with a gluino helps less or does not even help at 
all in reducing J7nlsp^ 2 - It can be seen from eq. ^ 
in Appendix A that, in the limit of a very strong de- 
generacy, (rag/mNLSP — 1) — ► 0, and a dominant gluino 
annihilation cross section, the resulting effective annihi- 
lation cross section is proportional to (16 + ,g)~ 2 , where 
.9 = 2 (6) is the number of the NLSP degrees of freedom 
for bino (wino) NLSP. Therefore we expect the effective 
cross-section for the wino to be smaller than for the bino 
and the wino abundance to be larger by ~ (11/9) 2 ps 1.5. 
Considering as well that the wino hadronic branching ra- 
tio is a factor tan -2 9w ~ 3.3 larger than for the bino, 
the stronger bound is explained. We see that a reheat- 
ing temperature of a few 10 9 GeV is also possible in the 
wino NLSP case with a 1% wino-gluino mass degener- 
acy, but due to the above effects such large reheating 
temperatures are possible only in the low mass region, 
for ma, < 200 GeV. Note that a light wino window, with 
masses around 100 GeV, is present also for the pure wino 
case, without coannihilations with the gluinos [47] . 

We also note that a 10% wino/gluino degeneracy ac- 
tually increases the NLSP relic density despite a larger 
annihilation cross section for gluinos, as it can be seen in 
Figure [9] as a little 'bump' in the predictions for ^nlsp^ 2 
in the wino NLSP case. This happens due to the presence 
of the 'weights' 7? in obeying Jo + j g = 1, and the 
fact that the real coannihilation cross-section involving 
a wino and a gluino in the initial state is negligible for 
larger squark masses. In this discussed in the ap- 

pendix A, the effective cross-section reaches a minimum 
when the ratio of the gluino over wino weights is equal 
to the ratio of the wino over gluino cross-sections. The 
increase in the abundance is at most 1 + <JqI Og , so it is 
negligible for the bino case, but visibile for the wino. 



C. Higgsino NLSP 

The case of higgsino NLSP in models of gauge media- 
tion is more involved. Although arranging for cancella- 
tions between various contributions to the soft sypersym- 
metry breaking Higgs mass parameter m 2 H , and hence /x 
parameter, is most welcome phenomenologically, it re- 
quires some fine tuning in the boundary conditions. Fur- 
thermore, for light neutralinos, the hadronic branching 
fraction is quite sensitive to the gaugino admixture and, 
e.g., 5% bino content can lower Bh even by an order of 
magnitude [47] . Moreover a non- vanishing gaugino frac- 
tion opens up also the channel of resonant annihilation 
via the pseudoscalar Higgs so that the higgsino number 
density can vary strongly around 2m^ ~ uia- For these 
reasons, we only indicate here the predictions for the 
higgsino NLSP in the conservative case of no resonant 
annihilation. With the gauginos twice heavier than the 
higgsino, we obtain the maximal reheating temperature 
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(2 — 3) x 10 7 GeV for higgsino masses between 100 and 
1000 GeV and 5% bino admixture. Larger reheating tem- 
peratures are surely possible if the annihilation proceeds 
on the resonance, at the cost of a fine-tuning between the 
neutralino and Higgs masses [4"T] . 

For the case of a higgsino NLSP, coannihilation with 
the gluinos is perhaps less natural since one could expect 
all the gauginos to be much heavier than the higgsinos. 
Also the equilibrium between higgsinos and gluino mainly 
proceeds through small Yukawa couplings and is some- 
what less effective than for gauginos. Finally, one needs 
to know the full spectrum of the higgsinos to account for 
the effective number of the degrees of freedom partici- 
pating in coannihilations. This can vary from 8, if the 
mass splittings among the higgsinos are not larger than 
with the gluino, to 2, if the mass gap between the lightest 
higgsino and the gluino is much smaller than the mass 
splittings among the higgsinos. The former case requires 
decoupling of heavy binos and winos: a situation with an 
extreme fine-tuning of the input parameters at the high 
scale of supersymmetry breaking and, as we have already 
seen in the wino NLSP large number of coanni- 

hilating states leads to larger Onlsp^ 2 ; hence, the latter 
case appears more plausible. 

We explored nevertheless the possibility of hig- 
gsino/gluino coannihilations within models of general 
gauge mediation with the messenger scale 10 15 GeV. We 
found an intermediate result between bino and wino neu- 
tralino: the annihilation cross-section lies between those 
for bino and wino NLSP, and in the limit of strong 
NLSP/gluino degeneracy one higgsino participates in 
coannihilations much more efficiently than the others. 
For a 300 GeV higgsino, we obtain the relic density of 
a few 10 -4 , while at the Higgs resonance it goes down to 
a few 10~ 5 . 
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FIG. 4: Mass spectra of sparticles in the example discussed 
in the text. The lowest neutralino line corresponds to = 
299.8 GeV, m N2 = 305.6 GeV and m Cl = 305.6 GeV, while 
Trig — 308.9 GeV. The Higgsino mass parameter is = 
1153 GeV and is irrelevant to the ensuing discussion. 



III. LHC DISCUSSION 

Taking into consideration the requirements of CP 
asymmetries, BBN, and gravitino dark matter, we have 
concluded that an interesting approach to viable leptoge- 
nesis leads to nearly degenerate gauginos and a spectrum 
of scalar superpartners somewhat higher in mass. In par- 
ticular, the gluino mass needs to be rather close in value 
to the NLSP mass such that coannihilation effects can 
sufficiently reduce the number density of the NLSP so as 
to not disrupt BBN when the NLSPs decay. This generic 
feature of the spectrum has important consequences for 
LHC discovery of supersymmetry. 



D. Gluino NLSP 



A. Gluino pair production 



Although gluino NLSP has a small relic density, with 
^nlsp^ 2 ranging from 6 • 10 -5 to 2 • 10~ 3 for gluino 
masses from 200 to 1000 GeV, successful primordial nu- 
cleosynthesis places very stringent limits on the presence 
of long-lived strongly interacting relic particles after the 
BBN [61 . If the NLSP is coloured, its lifetime should 
not exceed 300 sec and this is the origin of the constraint 
T R < 3 - 7 • 10 7 GeV for m g ranging from 200 to 1000 
GeV and the bino and the wino twice as heavy as the 
gluino. We also note here that the gluino NLSP relic 
density is smaller than the relic density of a stop with 
the same mass [62], hence the maximal allowed reheat- 
ing temperature in the stop NLSP case is lower than for 
the gluino NLSP. 



If the gauginos are much lighter than the scalars, then 
the largest supersymmetry production cross section at 
the LHC is gg. In the limit of a nearly exact degeneracy 
of the neutralino NLSP and gluino masses, there are no 
visible particles to trigger on from the production and 
decay of gauginos, since the final state is just an invisible 
neutralino and a very soft gluon or qq pair. For extreme 
degeneracy below a percent or two, the gluino decay is 
strongly suppressed and a displaced vertex at > 1 mm 
distances is possible, but again generally with no visible 
particle to trigger on, only one or two very soft jets. 

For illustration of the issues of detectability of super- 
symmetric particles at the LHC, we shall take a closer 
look at the model with mass spectrum shown in Fig- 
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channel 


branching fraction 


g -» N ig 


0.59 


g -> TVigg 


0.35 


9 -> N 2 g 


0.03 


g -> N 2 qq 


0.02 




0.41 


7V 2 -> iVi7 


0.31 


iV 2 -> iVil+l" 


0.08 



TABLE I: Branching fractions relevant for collider analysis. 



ure|4j This model can be realized within the framework 
of general gauge mediation with the messenger scale of 
10 15 GeV and tan/3 = 10. Bino/gluino degeneracy is 3% 
and the resulting maximal reheating temperature attain- 
able in this model is 3 • 10 8 GeV. From the branching 
fractions in table [I] we see that over 1 /3 of the gg events 
will be in the most advantageous channel gg — )• gNigNi, 
or in other words, two jets plus missing energy. 

For our example model point, the gluino mass is 
309 GeV and the leading order total cross-section at y/s — 
14TeV LHC is a{gg) ~ 255 ± 5pb. Thus, in a few ftT 1 
of data we expect quite a large number of events from gg. 
However, it will be difficult to trigger on these events and 
discern them above a large background. To see this, we 
note that in the rest frame of the gluino the energy of the 
gluon is fixed to be E g — (to? — m%)/(2m,g) = 9.0 GeV. 

We have conducted a MadGraph [63] Monte Carlo sim- 
ulation of the production of gluino pairs followed by de- 
cay into gluon plus neutralino. The results are given in 
Figure [5] where we have plotted the pr values of each jet 
at the parton level. Each event is a point in the {pT t , Pt 2 ) 
plane, where pr 1 is the higher px of the two jets. The 
gluon energy in the lab frame may increase or decrease 
depending on its relative decay direction to the boost di- 
rection. For this reason, the highest p? jet can be rather 
large - in this simulation of 1000 events, the highest pr 
obtained was nearly 60 GeV. 

Unfortunately, the visible energy and missing energy 
is not enough to trigger the events for saving at the 
LHC. There is too much background to open the trig- 
ger to events of this kinematic topology without sig- 
nificant prescaling that loses the signal. For example, 
for s/s = 7 TeV collisions, where the LHC is currently 
running, there are several triggers that are of poten- 
tial relevance to gluino gluino production followed by 
decays to soft things [53]; however, each fails to cap- 
ture a significant number of signal events. There is the 
"single jet" trigger requiring pr > 110 GeV, the "dijet 
trigger" which requires the average pt > 70 GeV, the 
"sum pt" trigger which requires p™ m (jets) > 200 GeV, 
and the non-prescaled "missing Et" trigger requiring 
MET > 60 GeV. These are only for the 7 TeV collider, 
and each of these numbers will be approximately doubled 
for the 14 TeV collider. Thus, none of these triggers will 



> 
a 




20 30 40 

Pti (GeV) 



FIG. 5: Scatter plot of pr values for the two gluons in 
gg ffA^iffiVi for 1000 events with m- g = 308 GeV and 
ulnx = 300 GeV. This simulation is for a pp collider with 
y/s = 14 TeV center of mass energy. 



efficiently record these events, and we seek a better path 
to discovery. 

When gaugino degeneracy is present, we have seen that 
the final states in the process are too soft to trigger on. 
Therefore, we need another process that can generate 
much higher pt jets or leptons. We remark that there is 
also the prospect of detecting gluino pair production via 
tagging from an initial state radiation (ISR) jet. This 
is a technique that has been advocated for many new 
physics scenarios that have no substantial visible energy 
from the final state [65] • ISR jets tend to be soft, and 
backgrounds are determined to large degree by how well 
the detector is understood and how well fake rates and 
jet energy measurements can be controlled. We do not 
pursue this approach here, but merely remark that it 
could be a useful signal for discovery, or even confirming 
a model if discovery is made through another channel. It 
would be especially important in the case of very high 
degeneracy such that the gluino decay has a displaced 
vertex, as discussed earlier in the section. 



B. Squark-gluino monojet signature 

We wish to determine if there might yet be another 
signature of value. Upon inspection of the spectrum it is 
evident that squarks are too heavy to pair produce effi- 
ciently, and certainly the same goes for the more weakly 
coupled sleptons. However, the squarks are often not too 
heavy for the promising signature of a single squark being 
produced in association with a (much lighter) gluino. 

An important contributing factor to the viability of 
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FIG. 6: Scatter plot of 1000 events at y/s = 14 TeV pp collider 
from qg — > qgg, where rrig = 300 GeV and = 1.5 TeV. The 
meaning of q is all first and second generation squarks: ul, 
ur, cLl, cLr, §l, sr, cl, cr. Each event is characterized by 
the absolute value of the pseudorapidity r\ of the quark jet 
from q — > qg decays, and the pr of that jet. Given that the 
gluino is nearly degenerate with all other gauginos including 
the NLSP, it acts like a source of missing energy. Thus, the 
Pt of the quark jet is approximately the missing energy of the 
event also. 



this signature is the high parton luminosities for gq ini- 
tial states that produce gq at tree-level. Furthermore, 
the kinematics of these events are favorable where the 
large mass difference between the squark and gluino gen- 
erates a large px jet from q — > gq decays. There is a 
also a large amount of missing energy recoiling against 
this jet. The signature is that of monojet plus missing 
energy. This has been studied recently within the con- 
text of the Tevatron for a wide class of models [55], as 
well as studies dedicated to the LHC, as we shall discuss 
below. Gluino pair production with ISR jet discussed in 
the previous section can also contribute to this signature, 
although the jet px is typically much softer. We shall ig- 
nore that additional small contribution to the high px 
signal defined in the following analysis. 

The model that we simulate is a slightly simplified 
version of our spectrum that we considered above. We 
consider here m g = 300 GeV and take the light squark 
masses to be 1.5 TeV in order to illustrate the signature 
with a reasonable spectrum. Generically, and certainly in 
this case, over 90% of light squark decays are into q —> qg, 
so single jet plus missing energy (i.e., soft-invisible decays 
of gluinos) becomes the most important signature to con- 
sider. The total leading order production cross-section 
at the 14 TeV pp LHC collider is 3.4 pb. Although this 
is nearly two orders of magnitude below the gg cross- 
section, the picobarn rate is high enough to record many 
thousands of events in the course of a few inverse fem- 



tobarns of integrated luminosity. Thus, it is a promising 
signature. In Figure [6] 1000 events are simulated, where 
we give the pseudorapidity (77) and the px values for each 
simulated event. To reduce backgrounds and increase re- 
liability of the analysis, it is often required that the jet 
be central (\n\ < 2) and have large px {px > 200 GeV). 
Those two requirements still leave the vast majority of 
events available for analysis. 

Figures [7] and [8] plot the total squark plus gluino pro- 
duction cross-section at a 7 TeV and 14 TeV center of 
mass energy pp collider as a function of the squark mass 
(first two generations) for various values of the gluino 
mass. When the squark mass is much greater than the 
gluino mass, the resulting jet px from squark decays to 
gluino plus quark are very high, and we get a strong single 
jet plus missing energy signature with very high trigger 
efficiency. 

The question of what the background is for the sin- 
gle jet plus missing energy is notoriously subtle. Our 
purposes here are to describe the basic features of the 
background, and give an estimate of expectations. Only 
a full detector simulation after careful engagement with 
the LHC data can ultimately determine what precise sen- 
sitivity levels can be reached. 

Nevertheless, we can compare our signal to the back- 
ground after cuts advocated in table 2 of [57] • First, 
the majority of our signal will be one jet plus missing 
energy. Multiple jets will arise from higher order correc- 
tions, which increases the signal; however, we do not in- 
clude these, thereby losing out on small additional signal, 
but also not being affected adversely by the "Number of 
jets < 3" cut implemented by [67 . Nevertheless, we know 
that that cut is very important in reducing tt and QCD 
backgrounds, and we can assume with overall impunity 
to the signal that it has been applied to the background. 
We can also apply E T USS > 400 GeV, which, given that 
we are working to leading order, automatically also im- 
plies the simultaneously required px(jet) > 350 GeV. It 
is automatic because at our leading order computation 
Px(jet) — E T USS to a good approximation. We can also 
require that the pseudo-rapidity of our jet is less than 
1.7 as |67| requires. The remaining two azimuthal angle 
cuts in table 2 of [67] have no consequence to us because 
they are automatically satisfied in our approximation. 
Furthermore, they have little effect on the larger back- 
grounds anyway. 

The combination of E^' 183 > 400 GeV and n jet < 1.7 
tends to reduce our signal by only ~ 30% which is within 
the uncertainty of QCD corrections and other uncertain- 
ties of the analysis. For the background, if we combine 
all these cuts we get approximately less than 50 events 
per 100 pb -1 of data. In other words, the background is 
about 0.5 pb. This estimate is also consistent with the 
results of refs. [HS1 EH] ■ Furthermore, these backgrounds 
can be measured well since at high missing Ex they are 
dominated by Z(yv) + j which can be normalized to the 
Z(l + l~) + j rate. Thus, assuming accurate computations 
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FIG. 7: Squark gluino production cross-section at \fs = 7TeV 
pp collider as a function of squark mass for various values of 
the gluino mass. The squark masses refer to first two gener- 
ation squarks. 



and measurements can be done for the background, and 
uncertainties can dip below 20%, we estimate reaching 
sensitivities down to approximately the 0.1 pb level for 
the signal. In 1 fb _1 of data, for example, we would have 
about 100 signal events at that cross-section, and results 
would be very near threshold for detectability. Thus, we 
can tentatively conclude that the 0.1 fb cross-section line 
is approximately the threshold for discovery of BSM sin- 
glet jet plus missing energy signature at 14TeV LHC. 
From that determination, everything above the dashed 
1 pb line in Figure [8] may be detectable at the LHC with 
over lfb -1 of data. Generically this leads to discovery 
sensitivity of gluinos less than ~ TeV with squarks above 
even a few TeV. 



C. Tevatron Remarks 

The parameter space we have been considering is where 
the gauginos, including gluino, are nearly degenerate. We 
should remark on what bounds there might be at the 
Tevatron for this scenario, if any. The limits on the gluino 
mass from CDF and DO experiments are usually quoted 
to be above ~ 300 GeV [70]; however, those limits assume 
that the gluino decays into an LSP of much lower mass 
such that the visible energy in the event is very high and 
can be triggered on. That is not the case in our model. 
The true capability of a limit should be much lower. 

When the gluino mass is degenerate with the LSP to 
within 10 GeV, the Tevatron has the potential to exclude 
its existence up to 160 GeV when the lumonisity exceeds 
2fb~ [7T]. The signal is gluino pair production in as- 



FIG. 8: Squark gluino production cross-section at y/s = 
14 TeV pp collider as a function of squark mass for various 
values of the gluino mass. The squark masses refer to first 
two generation squarks. The dashed line of a = 0.1 pb is the 
estimated lower limit of the cross-section that the LHC would 
be sensitive to when seeking a BSM contribution after more 
than lfb -1 of data is accumulated. 

sociation with a jet, leading to a monojet plus missing 
energy signature. Signal to background requirements are 
unlikely to allow improvement beyond 160 GeV even with 
greater luminosity [71]. It has been suggested [72] that 
the photon plus missing energy signature could be more 
probing than the monojet plus missing energy signature 
when luminosity exceeds lfb -1 . Although the signal is 
lower, the signal to background is higher for any given 
gluino mass, and with sufficient luminosity the higher 
signal to background signature always wins in sensitiv- 
ity. The probing sensitivity for gluinos nearly degener- 
ate with LSP in the photon plus missing energy chan- 
nel with greater than lfb -1 has been estimated to be 
175 GeV [75] ■ Further discrimination from background 
may be possible in some cases if gluino decays frequently 
produce leptons, albeit very soft ones, through interme- 
diate chargino or neutralino decays to the LSP. We are 
not aware of any Tevatron experimental paper that has 
reported an analysis for these scenarios, so the numbers 
above are suggestive of what could be done, and not what 
has been achieved yet. 



IV. SUMMARY & CONCLUSIONS 

In the case of a gravitino LSP and DM, we have ex- 
plored the parameter space of neutralino NLSP with a 
nearly degenerate gluino and we found that the BBN 
constraints are strongly relaxed in the bino case, so that 
a reheating temperature sufficiently high for thermal lep- 
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togenesis and the right gravitino abundance are obtained 
just with a mass degeneracy of the order of a few percent 
between bino and gluino. So in general a compressed 
gaugino spectrum helps in reconciling leptogenesis with 
gravitino DM and neutralino NLSP. On the other hand, 
for the wino and higgsino cases, the coannihilation with 
gluinos does not improve much the situation: for the 
wino case, it extends a bit the light wino window already 
found in [37] for a pure wino, while for the higgsino it 
does not beat the abundance suppression that can be ob- 
tained through the resonant annihilation. We found also 
that these degenerate gaugino scenarios can be embed- 
ded in general gauge mediation with a moderate tuning of 
the parameters, but since our results depend mostly only 
on the gaugino masses they are not restricted to these 
specific case and can be extended to any supersymme- 
try breaking scenario that allows for degenerate gaugino 
masses. 

In all these cases though, the gluinos are not only 
nearly degenerate with the NLSP, but also quite light 
< 300 GeV and are therefore copiously produced at the 
LHC. The fact that gravitino DM and thermal leptoge- 
nesis together give an upper bound on the gluino mass 
around the TeV scale was well-known [73] , but in our case 
in order to suppress the NLSP number density we need 
also a light NLSP and gluino. Unfortunately, since the 
degeneracy between them is small, even if the production 
cross-section is large, the main signal of a highly energetic 
track is missing and it is not so simple to disentangle the 
scenario from the QCD background. We suggested a cou- 
ple of final states and discuss if they could pass most of 
the LHC trigger cuts and some possible strategies for de- 
tection. The most promising channel is the gluino squark 
associated production, that gives a highly energetic jet in 
the final state and could be accessible even with lfb -1 
of data during the early phase of running of the LHC. 
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FIG. 9: NLSP relic density calculated in the presence of gluino 
coannihilations as a function of the degeneracy parameter 
Wg/mNLSP — 1 for rriNLSP = 300 GeV. Solid (dashed) lines 
correspond to bino (wino) NLSP. Thick (thin) lines show the 
results with (without) nonperturbative contributions. 



Appendix A: gluino coannihilations with 
Sommerfeld enhancement 



Here we discuss how to include nonperturbative ef- 
fects in the gluino annihilation, often referred to as Som- 
merfeld enhancement. Our calculation closely follows 

[31 Eg. 

With the NLSP of mass m, the freeze-out temperature 
is given by: 



x F 



In 



c(c + 2)i 



1 

27r 3 ^Jg*XF 



mMp(av) 



(3) 



where x = m/T, c is the coefficient in the relation 
Y = (1 + c)Y~ cq (which we set to 1/2), M P is the Planck 
mass (not reduced), g e g is the effective number of the 
degrees of freedom decoupling at freeze-out and (av) is 
the average containing the annihilation cross-section a. 
The average is: 



1 



(av) 



s 3/2«2 



T 



ds (4) 



with (3 = yjl - 4m 2 /s. 

We assume that the annihilation cross-section is the 
sum of the NLSP annihilation cross-section and the 
gluino coannihilation cross section: 



(5) 



where 70 = go/g c s, j g = 1 - 70, 

ffoff =50 + 16(1 + 8f' 2 t 



-5x 



(6) 
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S = TTig/mNLSP ~~ 1 arL( i 9o is the number of the degrees 
of freedom of the NLSP and equals 2 (6) for bino (wino) 
NLSP. This approximation is justified in the limit of 
heavy squarks, necessary to satisfy the Higgs mass bound 
with light gluinos. In this case in fact the coannihilation 
channel x + .9 ~~ ^ <7<7 is strongly suppressed by the large 
squark masses and can be neglected. Note on the other 
hand that the processes in which a NLSP scatters inelas- 
tically off a SM particle to become a gluino are efficient 
down to the freeze-out temperature, as the suppression 
of the cross-section due to heavy squarks is compensated 
by a large abundance of the light SM states. 

The neutralino annihilation cross-section includes 
many channels. For the gaugino neutralino case, an im- 
portant annihilation channel at low neutralino masses is 
the one into leptons via intermediate sleptons, which are 
in our case much lighter than the squarks. This annihi- 
lation rate is given by 



MA l/ r\ 2 
64-7T cos 2 6wP\s 



* + (1 + 2A) 2 - /?? 



In 



1 + 2A - ^ ' 
1 + 2A + /3* 



2A 



lz3 ' 

2(1 + 2A) 



(7) 



where A = (to? — m|)/s and 



A L = N jl ±N j2 cot9 w iQYi L jv u (8) 
A R = -2N n . (9) 



°{gg -> qq) 



160 g s 



(3-/?|)(3-/? 2 ) 



(11) 



with /3 a = \/T- 4m 2 /s. The contribution of the cross- 
section with intermediate squarks is negligible in this 
case. 

We see from the expression of the effective cross-section 
eq. ([5]), that the effective numbers of degrees of freedom 
7 change strongly with the temperature: defining as a 
variable the ratio z = 7^ /70 it is easy to see that 



a = a a 



1 + ^Vgg/vo 
(1 + z) 2 



(12) 



and one can show that the fraction on the r.h.s. has 
minimum as a function of z at z — o^jo-gg and that the 
minimal value is given as 



1 a ¥g 

1+ Z <Jgg + CT 



(13) 



So it is easy to see that this reduction of the cross-section 
is usually negligible for the case of bino neutralino, since 
Co <C <Jgg, but not for the wino. 

The Sommcrfcld enhancement is accounted for by mul- 
tiplying each of the cross-sections (10 1 and (111 by the 
factor 



Ei = 



1 — exp 



(_CiTTO JL \ 

V Pt J 



(14) 



For larger masses and especially for the wino case, also 
the channel into W + W~ becomes important and for the 
mixed case the resonant annihilation through the Higgs 
s-channel. In general in the case of wino neutralino also 
the coannihilation with the other charged winos is im- 
portant and the cross-section is not very much smaller 
than the gluino annihilation cross-section. Since we con- 
sider spectra for which the sfermions are generally heavier 
than the superpartner fermions, we use an approximation 
Co ~ P/s ( (To ~ l/ s /3) f° r bino (wino) NLSP, for which 
the p-wave (s-wave) annihilation dominates; the numer- 
ical coefficient is chosen so that it gives a correct relic 
density in the absence of coannihilations. In the case of 
the wino NLSP, one needs, in principle, to consider coan- 
nihilations between the wino states, but these states are 
typically extremely degenerate in mass, so uo should be 
considered an effective cross section with a common 70 
factor for all wino states in 

The gluino annihilation cross-section reads Ggg = 
cr(gg — > gg) + J2q a (99 ~ * where leading contribu- 
tions, mediated by the gluons, are: 



v{gg -> gg) 



3naj 
16pfs 



In 



1 + lg 
1-/35 



(21 



6$ 



3/?f)- 



33P g + 1701] 



(10) 



where Ci — 3/2(1/2) for the gg (qq) final state. The 
strong coupling constant is evaluated at scales Pgirig. 
This corresponds to averaging over the annihilating 
states, which produces a smaller correction than sum- 
ming up all the contributions. We also neglect the possi- 
bility of gluinos forming bound states, which would fur- 
ther enhance gluino annihiliations (often quite dramat- 
ically). Therefore, our assumption about including the 
nonperturbative effects in the gluino annihilations en- 
sures that the NLSP relic density calculated here is at 
worst an upper bound, so the claim that a given model 
satisfies the BBN bound is robust [52] . 

We can now calculate the final NLSP abundance from: 



^NLSP 



r 



dx 



(15) 



which can be used to evaluate the standard cosmological 
parameter 



^nlsp^ — 



2 ms 0^NLSP 



(16) 



where the present entropy density is sq = 2889.2 cm 3 
and the critical density is p c = 1.0539 • 10~ 5 GeVcm . 
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